Equation of state for pure SU(3) gauge theory on anisotropic lattices* 
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We present results for the equation of state for pure SU(3) gauge theory obtained on anisotropic lattices with 
the anisotropy £ = a a /a t = 2. The pressure and energy density are calculated on Nt/£ = 4, 5 and 6 lattices with 
the integral method. They are found to satisfy the leading 1/Nt scaling from our coarsest lattice N t /^ — 4. This 
enables us to carry out well controlled continuum extrapolations. We find that the pressure and energy density 
agree with those obtained using the isotropic plaquette action, but have smaller and more reliable errors. 



1. Introduction 

The study of lattice QCD at finite tempera- 
tures is important to clarify the dynamics of the 
quark gluon plasma [Q. In order to extract pre- 
dictions for the real world from lattice data, we 
have to extrapolate them to the continuum limit. 
However the continuum limit of the equation of 
state (EOS) has not been obtained in full QCD 
because of large lattice spacing dependence [|| Q . 

We have proposed using anisotropic lattices to 
solve this problem. As a first test, we have applied 
them for calculations of the EOS in pure SU(3) 
gauge theory and studied their efficiency 0. We 
report a summary of our results in this article. 

2. Simulation 

We employ the plaquette action given by 

s = P(r £ (! - P » £ (1 - 

so . 

n,ij n,i 

where £o is the bare anisotropy and (3 — Q/g$. 
We use the values of the renormalization fac- 
tor ?7(/3,C) = £/£o(/3,£) obtained by Klassen ||, 
where £ = a s /a t . 
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Our study in the Stephan-Boltzmann limit 
shows that £ = 2 is the optimal choice to re- 
duce discretization errors in the EOS §. There- 
fore we use £ = 2 anisotropic lattices with sizes 
JVf x N t = 16 3 x 8, 20 3 x 10 and 24 3 x 12. For 
N t — 8, we make additional runs on 12 3 x 8 
and 24 3 x 8 lattices to examine finite size ef- 
fects. The zero-temperature runs are made on 
N 3 x £N S lattices with £ = 2. After thermal- 
ization, we perform 20,000 to 100,000 iterations 
on finite-temperature lattices and 5,000 to 25,000 
iterations on zero-temperature lattices. 

3. Scale and critical temperature 

We determine the physical scale of our lattices 
from the string tension. The string tension a is 
extracted from the static potential V(R) at zero 
temperature assuming a form H], 



V(R) = V a +aR-e- + l 
R 



(2) 



where [1/R] is the lattice Coulomb term from one 
gluon exchange. Then we fit the string tension in 
(3 using an ansatz proposed by Allton M, 



a s ^ = f(J3) (1 + c 2 a({3) 2 + c 4 a(/3) 4 )/c , 



(3) 




where f((3) is the two-loop scaling function of 
SU(3) gauge theory and a[fi) = /(/3)//(/3 = 6.0). 

We define the critical coupling f} c {N t , N s ) from 
the peak location of the susceptibility for a 
Z(3)-rotated Polyakov loop. We extrapolate 
(3 c {N tl N s ) to the thermodynamic limit assuming 
a finite-size scaling law, 



(i c {N u N s ) = D (N t ,oo) -h(N t /(SN B ))- 



(4) 



From f3 c on anisotropic 12 3 x 8, 16 3 x 8 and 24 3 x 8 
lattices with £ = 2, we find h = 0.031(16) for 
Nt/€ = 4. We adopt this value for all N t . The 
critical temperature in units of the string tension 
is given by T c /^ = ^/N t a s ^ (0 c (N t , oo)). 

We extrapolate results for F = T c j ' y/a to the 
continuum limit with the leading scaling form, 



F\ = F\ 

\Nt I continuum 



c F /Nl 



(5) 



In the continuum limit, we obtain T c j \fa = 
0.635(10) from the £ = 2 plaquette action. 

Fig. [l] shows our results of T c /^/a. In Fig. |l|, 
we also plot the results obtained on isotropic lat- 
tices using the plaquette action [|| and the RG- 
improved action jl(J[ll). Our value of T c /^/a in 
the continuum limit is consistent with that from 
the £ = 1 plaquette action within the error of 
about 2%, but slightly smaller than that from the 
£ = 1 RG-improved action. This may be caused 
by the difference in the details of the potential fit. 



4. Equation of state 

We calculate the pressure p and energy density 
e using the integral method p2l, 
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The beta function d/3/da s \^ =2 is determined 
through the string tension parametrized by 
Eq. (||). We extrapolate the EOS to the contin- 
uum limit with Eq. (|^). 

Figs. U and |^ show the pressure and energy den- 
sity at T/T c = 2.5, as a function of (£/N t ) 2 (filled 
circles). For comparison, results from isotropic 
lattices using the plaquette action ]13|] (open 
circles) and the RG-improved action Ww (open 
squares) are also plotted. 

The advantage of anisotropic lattices is clearly 
seen in Figs. g an d [| On the coarsest lattice 
2Vf/£ = 4, the cutoff errors at £ = 2 are much 
smaller than those at £ = 1 with the same plaque- 
tte action. Comparing the computational cost to 
achieve comparable systematic and statistical er- 
rors on isotropic and £ = 2 anisotropic lattices, we 
find that the anisotropic calculation has a factor 
of approximately 5 gain. Furthermore, since the 
£ = 2 data satisfy the 1 /iV t 2 scaling from N t /£ = 4 
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Figure 3. e/T 4 at T = 2.5T C . 
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Figure 4. p/T 4 and e/T 4 as functions of T/T c . 



(the right-most point), we can reliably perform an 
extrapolation to the continuum limit using three 
data points. Results at other temperatures are 
similar. 

Our results for the EOS in the continuum limit 
are shown in Fig. |J. The EOS from £ = 2 is 
consistent with that from the £ = 1 plaquette 
action within the error of about 5%. Notice that 
our £ = 2 results have smaller and more reliable 
errors of 2-3%. On the other hand, the EOS from 
the £ = 1 RG-improved action is larger by 7-10% 
(about 2a). than that from our £ = 2 results 
in the continuum limit. A possible origin of this 
discrepancy is the use of the iV t /£ = 4 data of the 
RG-improved action, which show a large (about 
20%) deviation from the continuum limit. 

5. Conclusion 

We have studied the continuum limit of the 
EOS in pure SU(3) gauge theory on £ = 2 aniso- 
tropic lattices, using the plaquette action. Aniso- 
tropic lattices are shown to be efficient in reduc- 
ing the cutoff dependence of the EOS. As a re- 
sult, we can perform continuum extrapolations 
from our coarsest lattice AT t /£ = 4. The EOS in 
the continuum limit agrees with that obtained on 
isotropic lattices using the same action, but has 
much smaller and better controlled errors. 

Anisotropic lattices may help us extract the 
continuum limit of the EOS in full QCD. 
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